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Abstract 

Under the NASA Fundamental Aeronautics Program the Su- 
personics Project is working to overcome the obstacles to super- 
sonic commercial flight. The proposed vehicles are long slim 
body aircraft with pronounced aero-servo-elastic modes. These 
modes can potentially couple with propulsion system dynamics; 
leading to performance challenges such as aircraft ride quality 
and stability. Other disturbances upstream of the engine gener- 
ated from atmospheric wind gusts, angle of attack, and yaw can 
have similar effects. In addition, for optimal propulsion system 
performance, normal inlet-engine operations are required to be 
closer to compressor stall and inlet unstart. To study these phe- 
nomena an integrated model is needed that includes both air- 
frame structural dynamics as well as the propulsion system 
dynamics. This paper covers the propulsion system component 
volume dynamics modeling of a turbojet engine that will be 
used for an integrated vehicle Aero-Propulso-Servo-Elastic 
model and for propulsion efficiency studies. 

Nomenclature 

A area, m 2 

c p specific heat at constant pressure, J/(kg*K) 

c v specific heat at constant volume, J/(kg*K 

F thrust, N 

g gravitational constant, 1 (kg*m)/(N*sec 2 ) 

h static enthalpy, J/kg 

/ moment of inertial, (N*m)/sec 2 

J mechanical equivalent of heat, 1 (N*m)/J 

K a , K b combustor coefficients, unit less 
K c combustor coefficients, (N 2 *sec 2 )/(kg 2 *m 4 *K) 

K b bleed flow coefficient, (kg*K 1/2 )/(N*sec) 

K noz nozzle parameter, (kg*m 2 *K 1/2 )/(N*sec) 

/ length, m 

M mach number 

M air molecular weight of air, 0.02897 (kg/mol) 

Mjuei molecular weight of fuel (JP-4), 0.139 (kg/mol) 

M ca molecular weight of cooling air, 0.02897 (kg/mol) 

N rotational speed, rpm 

N c corrected speed ratio 

P pressure, N/m 2 

P r pressure ratio 


R universal gas constant, 287 (N*m)/(kg*K) 

r rotor mean radius, m 

r T rotor tip radius, m 

T temperature, K 

U velocity, m/sec 

U ab volumetric flow rate of combustor air, m/sec 
V volume, m 3 

W mass flow rate, kg/sec 

Wcmf corrected mass flow rate, kg/sec 

Wf externally acted upon fuel flow, kg/sec 

Greek 

y ratio of specific heats, y cp = 1 .4, y cb = y tb = y ab = 1 .3 1 

r\ efficiency 

0 molecular vibration energy constant associated with a 
gas mixture, K 

p weight density, kg/m 3 

t delay time 

T t delay time associated with fuel transport and mixing 

§ fuel to air ratio 

Subscripts 

a variable associated with combustor acoustics 

ab variable associated with the afterburner 
amb variable associated with ambient conditions 
b variable associated with compressor bleed 

c variable associated with stage characteristics 

ca variable associated with cooling air 

cb variable associated with the combustor 

cp variable associated with the compressor 

d variable associated with designed value 

e variable associated with engine exit condition 

/ variable associated with fuel 

fl variable associated with flame dynamics 

1 index associated with nozzle volume element number 

j index associated with turbine stage number 

k number of compressor stages 

/ index associated with inlet volume element number 

m number of turbine stages 

n index associated with compressor stage number 

noz variable associated with the nozzle 
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o variable associated with engine entrance condition 

q number of inlet volume elements 

perf variable associated with a thermally perfect gas 
ref variable associated with reference value 
tb variable associated with the turbine 

s static condition 

sv static condition variable associated with stage volume 

tc total condition variable associated with stage charac- 
teristics 

tv total condition variable associated with stage volume 
z variable associated with nozzle 

Introduction 

In many aero-propulsion applications a relatively low fidelity 
engine model (i.e., without volume dynamics), which includes 
speed and temperature dynamics with the appropriate component 
performance characteristics is adequate to design control laws. 
The actuation system, such as fuel actuation, is rather slow and as 
a result the control bandwidth is reduced, which makes such 
lower fidelity simulations (refs. 1 and 2) appropriate for controls 
design. In supersonics, the slim body structure of the vehicle 
excites structural dynamics, which can introduce upstream flow 
field disturbances at higher frequencies (up to 50 Hz or more). 
The vehicle control surfaces are actuated to suppress some of 
these structural modes associated with flutter (refs. 3 and 4). 

Analysis needs to be carried out to understand how these aero- 
servo-elastic (ASE) excitations as well as other atmospheric 
disturbances impact the propulsion system in terms of thrust 
variations that can affect ride quality and vehicle stability. In 
addition, the coupling of the propulsion system to the vehicle 
ASE modes needs to be analyzed, and for that an integrated aero- 
propulso-servo-elastic system (APSE) simulation is needed. Due 
to the high frequencies of these modes the propulsion system 
simulation needs to incorporate ID volume dynamics based on 
conservation equations modeling. 

A propulsion system for a supersonic vehicle also differs from 
a conventional vehicle in that a supersonic vehicle utilizes a su- 
personic inlet with active controls for shock positioning, which 
can also excite higher frequency dynamics. In addition, how these 
upstream disturbances impact the inlet shock positioning, pres- 
sure recovery, as well as inlet distortion and its impact on thrust 
variations also needs to be analyzed. Also to achieve high per- 
formance, a supersonic vehicle would operate closer to compres- 
sor stall and inlet unstart. This necessitates modeling utilizing 
component volume dynamics in order to study the dynamic inter- 
action of these engine components. This paper covers the propul- 
sion system component volume dynamics modeling of a turbojet 
engine. The supersonic inlet dynamics and their effect on the rest 
of the propulsion system and the vehicle will be addressed in 
some future studies. 

The propulsion system simulation described in this paper uses 
the architecture of a J85-13 turbojet engine; although it can be 
adapted to simulate different engines due to its generic architec- 
ture. The development of this engine simulation in MATLAB 
SIMULINK® is partially based on a simulation developed years 


ago at NASA Glenn Research Center (GRC) (ref. 5), using an 
analog computer and also referred herein as the Seldner simula- 
tion. In the engine simulation described here, some of the geome- 
tries are approximated, and the component performance maps are 
generic. The Seldner simulation was verified in both steady-state 
and dynamic modes by comparing analytical results with experi- 
mental data obtained from tests performed at GRC with a J85-13 
engine. The simulation described in this paper has been verified 
by comparing it against simulation results provided in Seldner 
(ref. 5), and also from simulations that display expected engine 
responses based on experience. 

The engine model presented is based on stage-by-stage com- 
ponent volume dynamics. The simulation results presented are 
based on component level volume dynamics. The component 
level volume dynamics model of the engine is envisioned to be a 
stepping stone towards eventually obtaining a stage-by-stage 
model with all the appropriate stage performance maps that may 
be required to address the thrust variations issue. 

The paper is organized as follows. First, the engine dynamic 
model development is described in detail for each major compo- 
nent. This is followed by some simulation results to help verify 
the engine simulation, as well some preliminary thrust variation 
studies. Finally, future plans for upgrading the engine model, 
including plans for additional analysis and controls design are 
discussed, followed by concluding remarks. 

Engine Model 

The Seldner (ref. 5) simulation technique was applied to the 
J85-13 engine. The simulation developed herein is also applied to 
the J85-13 engine, with some assumed component performances 
and geometries. The J85-13, shown in figure 1 is an afterburning, 
turbojet engine with an eight-stage axial-flow compressor and a 
two stage turbine. Compressor discharge air that bypasses the 
combustion zone is used to cool combustor and turbine compo- 
nents. The engine also has variable inlet guide vanes, controlled 
interstage bleeds, and a variable-area exhaust nozzle. The model- 
ing approach described herein is generic. But the reason this 
engine was chosen is that some information and simulation re- 
sults about this engine are available, and also because this engine 
architecture is closer to an engine architecture that is expected to 
be used for a supersonic vehicle. 

The engine simulation incorporates component performances 
in the form of performance maps for the compressor and turbine 
along with component volume dynamics. The modeling ap- 
proach follows the basic approach discussed in reference 5 for 
the J85-13 engine simulation, developed utilizing an analog 
computer. Differences in the modeling approach are noted in the 
upcoming sections of this paper that cover component simula- 
tions. Some of the pertinent component performance maps and 
geometry properties of the engine were not recorded in refer- 
ence 5. Therefore, some assumptions are made in developing 
this simulation. The simulation is verified against steady state 
experimental results and some dynamic responses available in 
reference 5. The simulation described in reference 5 covers 
stage-by-stage control volume dynamics, while the actual 
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Exhaust nozzle 



Figure 1. — J85-13 engine. 


simulation here is incorporated with component level lumped 
volume dynamics, and therefore, there could be some differ- 
ences in dynamic responses. 

Compressor 

This section describes the compressor model based on the 
individual compressor stages. In reference 5 the compressor 
maps were expressed in terms of normalized coefficients of 
pressure and temperature versus flow. Because it is rare at 
present to find performance maps expressed in these terms, a 
map generation routine is used to derive overall compressor 
maps in a more conventional form consisting of the pressure 
ratio versus corrected mass flow rate and compressor effi- 
ciency versus pressure ratio. The compressor maps would 
normally differ for inlet guide vane position and bleed flow. 
Certain geometrical properties of the compressor like cross 
section area, length, mean radius, and tip radius are not avail- 
able and are estimated. 

For a stage-by-stage compressor model, the individual com- 
pressor stage characteristics are needed. This includes comput- 
ing the pressure and temperature rise across each stage, by 
representing each stage with a pair of pressure and temperature 
(efficiency) maps, followed by its appropriate stage volume. 
Based on such an approach, the compressor consists of inter- 
connected stages, stacked together to form the overall compres- 
sor model. An idealized compressor stage is shown in figure 2. 

The gas dynamics associated by applying continuity, mo- 
mentum, and energy to a compressor stage, are 


n-1 stage ■ 


■ n stage- 


■ n+1 stage 



tage -HI 


'W r 


tv,n . 


A n , V n , l„ 


n th stage 
Compressor/ \ 'stator 
element / \ rotor 


'W f 


c,n+l 


n th stage equivalent 
lump volume 


Figure 2. — Schematic of n th compressor stage. 
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Yep 


i. tv,n Tr 

dt Psv,n^n 


( T tc,nK,n - T tv,nK,n + 1 - T tv,nK,n ) ( 3 ) 


where total temperature and pressure conditions at the com- 
pressor first stage and inlet exit are equivalent (T tc>n=1 =T t> i= q 
and P tc>n =i=P t> i= q ), and mass flow at the compressor exit and 
combustor inlet are equivalent (W c>n =k+\=W c b). The equation of 
state for the n th stage volume, for k number of stages, is 


r 


P - 

1 tv,n 


i + y spu, 2 


V 


Yep " 1 


P sv,n^^tv,n 


( 4 ) 


J t Ps V ,n=^{Kn-K,n + l-Kn) 


( 1 ) 


— \y = (p ~p ) 

i. c,n j V* tc,n ± tv,n-f 

at L 


f y -1 A 

Yc P hj 2 


1 + - 




Y cp - 1 


( 2 ) 


The correction factors involving Mach number and gamma in 
eq. (2) and (4) are inserted to correct for the differences be- 
tween static and total conditions. If these corrections are not 
made, the steady-state accuracy of the simulation will not be 
impacted, but dynamically the error can be up to 15%. The 
stage Mach number is 
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M, 


n 

P 


sv, n 


w 

rr c,n 

^ri yjlcp RT s,n 


(5) 


where 


T -T - 

s,n tv,n 


WS 


2 p; 


A 2 c 

sv,n n p,n 


( 6 ) 


Bleed flow is extracted from some of the compressor 
stages. The bleed flow relation assumes choked flow condi- 
tions, were K b , is the portion of the bleed flow effective area. 


W bn =K h A i 


b^b,n 


If 

A tv, n 


( 7 ) 


The overall compressor performance maps are shown in 
figure 3 and figure 4. For stage-by-stage modeling, individual 
stage maps should be utilized instead. The corrected mass 
flow rate is shown in figure 3. If this is the map for the n th 
stage, the corrected parameters are computed as 


W, 


cmf,n 


w If 

Tr c,n y tv, n- 


■1 j R ref 


R tv,n - 1 / R ref 


( 8 ) 


Similarly, the corrected speed ratio used in figure 3 and 
figure 4 is computed as 

(9) 

7V d,n 

The corrected mass flow rate and corrected speed are used to find 
the pressure ratio in figure 3, and the pressure ratio and corrected 
speed are used to find the efficiency in figure 4. From the pres- 
sure ratio, P rm the stage total pressure can be computed as 

R tc,n ~ R r,n * R tv,n - 1 (10) 


From the efficiency ratio and the pressure ratio, the stage total 
temperature can be computed as 


T — T 

± tc,n ± tv,n-\ 


1 + 


pftcp 
1 r,n 


\cp 


The cross section area in eqs. (2), (5) and (6) is 

A n = 4nr n( r Tn ~ r „ ) 


( 11 ) 


( 12 ) 


The differences between the model presented here and the 
one from reference 5 are as follows. First, in reference 5, cor- 
rection factor approximations are given for the momentum and 
the state equations but are not utilized in the simulation. In- 
stead, in this development exact relationships are provided. 



0 10 20 30 40 50 

Corrected Mass Flowrate ( W cfm ), [Ibm/sec] 


Figure 3. — Overall compressor pressure ratio map. 



Figure 4. — Overall compressor efficiency map. 


Second, the compressor maps are presented here in terms of 
pressure ratios and efficiency, instead, of flow coefficients that 
were utilized in reference 5. Third, based on these changes, all 
the pertinent mathematical relations for the Mach number and 
the corrected map quantities are provided in eqs. (5) to (11), 
including the cross section area calculation which is provided in 
eq. (12). 

Equations (1) to (12) represent the stage modeling of the 
compressor. In the absence of individual stage maps, the root 
of the number of stages of the compressor pressure ratio and 
efficiency ratio (i.e., the parameter raised to the power of 1/n) 
could serve as a rough approximation. However, the accuracy 
of the compressor and turbine simulations results and stability 
could depend on the accuracy of the stage or component per- 
formance and geometries. Stacking the stages together and 
with the individual stage performance maps, the compressor 
model is complete. 

Combustor 

The combustor is lumped into a single equivalent one- 
dimensional volume. In this volume the fuel flow combustion 
dynamics are also added. 

The combustor geometry is chosen to represent a lean burn- 
ing combustor, which is more appropriate for future combus- 
tors. Typically, in a behavior type model of lean burning 
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combustors (ref. 6), the combustion dynamics can be modeled 
as a self excitation system of a first order transfer function 
(TF) representing the flame dynamics, with a second order 
undamped TF representing the acoustics. But the acoustics 
here are represented with a first order TF. The assumption is 
that the unsteady combustion typical to lean burning combus- 
tors will be mitigated by some control approach like the one 
described in (ref. 6). The total combustion time delay is the 
sum of the delays of fuel transport and mixing, the flame dy- 
namics, and the acoustics as 

*cb=*t+*fl+*a ( 13 ) 


tors are not used in the combustor simulation due to the low 
Mach number in the combustor, which makes a negligible 
difference. The c p value in eq. (18) is not constant in the upper 
temperature range of the combustor, and for more accurate 
calculations eq. (44) that is provided later can be used. 

In this development, h c was assumed to be a constant of 
42.8e6 [J/kg] for JP-4 fuel, at about 1000 K (ref. 7). The com- 
bustor efficiency map (ref. 5), is shown in figure 5, where 

A T t,cb = T t,cb ~ T tv,n=k ( 2 1 ) 

The pressure drop across the combustor is 


Given a combustor volume, V cb , and a volumetric flow rate 
through the combustor, U cb , the combustor time delay can be 
approximated as 


• 2 

A P t , cb = + K B T t , cb ) 

*tv,n=k 


( 22 ) 


x cb=Vcb/U cb (14) 

The total combustion time delay is assumed to be in the order 
of 5 msec, with the fuel transport and mixing time delay ac- 
counting for the greater part of this delay. The overall fuel 
flow combustion dynamics are modeled as follows 


W f 


Wf 


Ke~ x ‘ s 


(x^5 + l)(x a 5 + l) 


(15) 


where K is a proportional TF gain, set to one. 

The combustor gas dynamics for continuity, momentum, 
and energy are 

j t 06 ) 

= 4 ~ fe-t -6U -"U) (17) 

dt leb 


where K /t and K B are experimentally determined; K A is found 
from non-combustion flow tests; K B from combustion flow 
test. The proportionality constant, K c , is found by solving 
eq. (22) for K c , given the designed steady-state quantities in 
eq. (22) at different speeds. A table is provided in reference 5 
for these values. Typical values for AP Ucb range from 0.05 to 
0. 1 . In reference 5, K c has the same units as K A and K B , but K A 
and K b should be dimensionless. 

The differences between the combustor model presented 
here and the Seldner model, reference 5, are as follows. First, 
the time delays in eqs. (13) to (15) are chosen for lean burning 
combustors and the relation of the overall combustor time 
delay, eq. (14), is provided here. Second, the equation of state 
is given here, where the specific gas constant is formulated for 
the combustion mixture. Third, the enthalpy, h c , was used as a 
constant here of about 1000 K. In the Seldner model, h c was 
determined using stage characteristics with the independent 
variable T t>cb . Fourth, a correction is made to eq. (22) here that 
involves W cb computed from the combustor volume dynamics 
instead of the mass flow rate coming from the compressor. 

Equations (13) through (22) complete the combustor model. 


ji _ y cb 

dt t ’ cb ~P s ,cbKb 


T tv,n=k^cb + W f T t,cb~ T ty cb^c,j=\ 


+ 3c^w h 

C p,cb 


and the equation of state is 

^t,cb ~ Ps,cb-^cb^t,cb 


R cb ~ 


R 


(l air j- ue i 


(18) 


(19) 

( 20 ) 


where a fuel/air ratio, 4), of about 0.02 was used that represents 
lean type combustion. Unlike the compressor, correction fac- 



Figure 5. — Combustor efficiency representation. 
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Turbine 

In the Seldner model, reference 5, the entire turbine is mod- 
eled as a single lumped volume. Moreover, the turbine is 
combined with the non-combusting afterburner and the nozzle. 
This combined model was the result of having the turbine 
contributing to the energy conservation through the enthalpy 
change, the afterburner contributing to the momentum by 
dominating the total volume, and both the turbine and the 
afterburner sharing the continuity equation, with the nozzle 
governed by compressible flow and choked flow conditions. 
While the model in reference 5 may exhibit sufficient steady- 
state accuracy, it’s not anticipated that dynamically the fre- 
quencies of these components will be accurately represented. 
Therefore, the turbine and the subsequent afterburner models 
are different than those presented in reference 5. 

An idealized turbine stage would be similar to the compressor 
stage shown in figure 2. The gas dynamics associated by applying 
continuity, momentum, and energy to a turbine stage are 






—T = 
dt ^ 


P sv,i V i 


(TtcAj + T t,caW caJ - T^jW^ 1 ) (25) 


where total temperature and pressure conditions at the turbine 
first stage and combustor exit are equivalent ( T tc j =1 =T t>cb and 
Ptcj=l ~ P i, cb ), and mass flow at the turbine exit and afterburner 
inlet are equivalent ( W c j =m+X - W ah ). The equation of state 
for the / h stage volume is (for m-number of stages) 



12 13 

Corrected Mass Flowrate [Ibm/sec] 

Figure 6. — Overall turbine pressure ratio map. 



7 8 9 10 11 12 13 14 

Corrected Mass Flowrate [Ibm/sec] 

Figure 7. — Overall turbine efficiency map. 

The c p value in eq. (29) is not constant in the upper tempera- 
ture range of the turbine, and for more accurate calculations 
eq. (44) that is provided later can be used. 

Performance maps are incorporated in this development that 
are more readily available, such as the ones shown in figure 6 and 
figure 7 for the overall turbine pressure ratio and efficiency. 

The corrected mass flow rate and speed ratio are used in 
figure 6 and figure 7. If these are the maps for the / h stage, the 
corrected parameters would be computed as 


p »j= 7 - 


] r *~ l p sv n R tb Ttv, n (26) 
2 ) 


W cmf ,j 


V/ ^ T tv, j-1 ^ P ref 

PfvJ-1 / Pref 


where 


N c,j = T7" 
J Nj 


\ T d,j 


(1 - <\))(M air + M ca ) + §Mj uel 


From the pressure ratio, P rj , the stage total pressure can be 
computed as 

Ptcj = PtvJ- 1 / p r,j (32) 


psv,f A fJy t b R tb T s, 


From the efficiency ratio and the pressure ratio the stage total 
temperature can be computed as 


T . - T • 

1 S,J 1 tv,j 


2Psv,j A j c pj 


r„ , =- 


l ~ P r 
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The cross section area in eqs. (24), (28) and (29) is 

A j =4nr j (r Tj -r j ) (34) 

Equations (23) to (34) represent the stage modeling of the 
turbine. Stacking the stages together and with the individual 
stage performance maps, the turbine model is complete. 

Afterburner and Nozzle 

A combined combusting afterburner and exhaust nozzle 
model is presented here, with the assumption that a supersonic 
cruise vehicle would likely employ an afterburner. It is assumed 
here that the large afterburner acts as a large filter, attenuating 
high frequency upstream disturbances, such that a detailed 
multi-finite element volume nozzle model will not impact dy- 
namic thrust calculations. With this assumption, in the com- 
bined model the afterburner volume dominates. Also, unlike the 
Seldner model, a combusting afterburner is modeled here, 
which also allows for temperature changes across the after- 
burner. Like the Seldner model, the nozzle is considered as a 
variable area compressible flow passage capable of choking. 

Based on these assumptions, and by separating the after- 
burner from the turbine model as discussed in the turbine 
section, the gas dynamics associated by applying continuity, 
momentum, and energy to the afterburner and nozzle are 

s,ab =^(Wab + W fuel , ab ~W Z ) (35) 

dt V ab 


The second and third multiplicative terms in eq. (39) represent 
the well known compressible flow function (ref. 5). The pa- 
rameter K z is variable and proportional to the nozzle area. It 
represents the turbojet terminal impedance, and lumps the 
nozzle flow coefficient and nozzle area with other flow pa- 
rameters. This parameter, through the variable nozzle area, is 
scheduled as a function of engine speed to match the steady- 
state operating line. In the simulation, K z was experimentally 
determined to better match expected results. 

Equations (35) to (39) complete the afterburner and the 
nozzle model. 

Rotor Dynamics 

The steady-state performance of a turbojet engine matches 
the compressor with that of the turbine operating points. A 
mismatch in these components produces an unbalanced torque 
or acceleration, which is integrated through the dynamic rela- 
tions to seek a new steady-state match. The rotor dynamics are 
based on the changes of mass and enthalpy. Therefore, the 
change of rotor speed is a function of the energy differential 
between the work extracted by the turbine and the work done 
by the compressor a 

T< n = f— ] 4^[(Kj=Ab + w ca Ka - w ab h j=m ) 

dt v 7i J NI LX J J ’ (40) 

“ {^c,n=kK=k +W b h b - W c n=l h n=1 

The enthalpies in eq. (40) can be computed as follows. 


±W h= ^ab8_ 

dt l ab 


( P tv,j=m ~ P t,ab)\ 1 + 


Vab 

Yah ~ ^ Vab~l 


ab 


(36) 


-\ ab =-^(T tvJ=m W ab +Q ab -T Uab W z ) (37) 

dt Ps,ab V ab 


At some future point, the heat addition Q ab will be consisting 
of two terms, similar to the combustor model. One term will 
be due to the heat addition of the afterburner fuel, and the 
other term will be due to the heat addition of the combusted 
fuel mixture with the associated fuel system dynamics. The 
equation of state is 


^ cb c p,cb P t,cb 

(41) 

ev, 

II 

(42) 

h - r T 
,l n ° p,n ± tv,n 

(43) 


The specific heat at constant pressure, c p , is not constant and it 
varies as a function of temperature and gas mixture based on 
the following equation, for a thermally perfect gas. 


c p,x ( c p ) perf 


1 + 


Y perf 


Y perf 


®x 

\Jt,x J 


e v /r M 


{e x ' T ^ -l ) 2 


(44) 


7 

p f i . Tab 1 ^ Yab~l n 

F t,ab- l V + ~ 2 M ab I Ps 


, R 


where R ab can be calculated similarly to eq. (20). The nozzle, 
with its mass flow rate, becomes the terminal boundary of the 
turbojet as 


= 


K z P t,ab 


vs 


ab 



1 



p A 

± s,amb 

V ab 

\ 

1- 

p , 

± s,amb 

v P t,ab ) 

1 


[ P l,ab j 


hid 

V ab 


(39) 


The combustor and the turbine specific heat vary some in the 
order of 1000 K. But the specific heat will be different due to the 
addition of fuel. For a pure air mixture, 0 = 3056 K. 

The differences between the rotor model presented here and 
the Seldner model, reference 5, are as follows. First, the Seld- 
ner model does not take into account the cooling air coming in 
the turbine. Second, in the Seldner model the enthalpies are 
not defined, and third, that model assumes a constant specific 
heat. 
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Engine Thrust 

Thrust was not included in the Seldner model, but it is esti- 
mated in this model as follows. In general, the net thrust is the 
sum of the momentum thrust and the pressure thrust as 

W 

F ~ —^~(U e - U Q ) + (P e - P SAmb )A e (45) 

g 

and 

Yz 


At steady-state, the engine thrust is primarily due to the momen- 
tum thrust; the pressure difference in eq. (45) is about zero by 
scheduling K z (i.e., by adjusting the nozzle exit area) as a function 
of the speed, N. However, in the dynamic sense, for thrust calcu- 
lations, the pressure thrust in eq. (45) could be significant. 


F e - Pt,ab 


1 + - 


!-Yz 


hz ( 1 + Yz) 


Engine Simulation Results 

In this section, selected results are presented for the whole 
engine simulation running together at sea level static flight 
conditions. These are for the component level lumped volume 
dynamics engine model. All it takes to convert the compressor 
and turbine models presented earlier into component level 
lumped volumes is to drop the subscripts that represent the 
stage-by-stage model and also use the performance maps for 
the entire component. Some of the initial objectives are to 
analyze thrust variations and controls designs using the com- 
ponent level lumped volume dynamic model, and later to 
address whether a more detail model is needed that employs 
stage-by-stage volume dynamics. 

Table 1 shows some steady state comparisons at sea level 
conditions between the engine, the Seldner simulation, and the 
model developed in this paper. The simulation results closely 
depict the engine steady state performance, within a few per- 
cent errors. One point worth noting here is that the engine 
model in this paper employs component maps and geometries 
that are not the same as the ones used in the Seldner model. 
Therefore, direct comparison of simulation results is difficult 
and can only be done in a relative sense. But even that would 
be challenging because the engine operating conditions that 
were used to produce these results, for Seldner model, were 
not fully documented. 

TABLE 1.— COMPARISON OF STEADY STATE 


RESULTS AT SEA LEVEL 



Rotor 

Speed 

Compressor 

Massrate 

Turbine 

Pressure 

Turbine 

Temp. 

Experimental 

100% 

19.90 [kg/sec] 

246,200 [Pa] 

992 [K] 

Seldner Model 

100% 

20.00 [kg/sec] 

264,100 [Pa] 

1021 [K] 

This Model 

104% 

20.99 [kg/sec] 

248,320 [Pa] 

955 [K] 

Error (Exp. Vs. 
this model) 

4.00% 

5.48% 

0.86% 

3.72% 


Figure 8 and figure 9 show comparisons of frequency re- 
sponses from simulation results obtained from the model in this 
paper and results from the experimental engine published in 
Seldner’ s report. These frequency responses show that the magni- 
tudes are about the same, with some variation in the frequency 
range of 1 to 40 Hz. Also, the dominant frequencies in these 
responses are nearly the same, in the range of about 35 to 50 Hz. 
Simulation results from the Seldner model, not shown here, more 
closely approximate those of the experimental engine. The reason 
is that the Seldner model matches the experimental engine in 
terms of geometries and component performance. 

Figure 10, figure 11, and figure 12 show engine step re- 
sponses due to the combustor fuel mass addition, at the same 
steady state operating point depicted in table 1 . These responses 
follow trends that would be expected from this engine. 
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Figure 8. — Comparison of engine simulation and experimental 
results for compressor pressure due to fuel response. 



Figure 9. — Comparison of engine simulation and experimental 
results for combustor pressure due to fuel response. 
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Figure 12. — Step response of turbine temperature to fuel. 



Figure 15. — Step response of compressor total 
temperature to engine input total temperature. 
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Figure 16. — Step response of turbine total 
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Figure 13. — Step response of combustor 
pressure to engine input total pressure. 


temperature to engine input total temperature. 




Figure 14. — Step response of turbine total 
pressure to engine input total pressure. 

Figure 13 and figure 14 show engine total pressure step re- 
sponses due to engine input total pressure at the compressor 
face, again at the same operating condition as shown in table 1 . 
The input pressure is stepped by 250 Pa, then back to nominal, 
followed by a -250 Pa step, and back to nominal. 

Figure 15 and figure 16 show engine total temperature step 
responses due to 1 K engine input total temperature steps at 
the compressor face. These responses also follow expected 
trends, with the turbine temperature step response being a bit 
more complicated. It has been observed that the turbine tem- 
perature response is rather sensitive to the accuracy of the 
turbine pressure ratio map. Thus, the turbine temperature 
transient response may change somewhat in the future with the 
substitution of more precise maps. 


Figure 17. — Frequency response of engine thrust due to engine 
input total pressure with different afterburner lengths. 

Figure 17 shows the TF of the thrust response due to the en- 
gine input total pressure for two cases; an afterburner length of 
0.75 m and an afterburner length of 1.5 m (nominal). As seen in 
this figure, the effect of halving the afterburner length (halving 
its volume) is to shift the frequency response to the right, ap- 
proximately double the frequency. This indicates that for this 
engine architecture, the afterburner dominates the frequency 
response due to its large volume in comparison with the rest of 
the engine components. These are preliminary results. But if 
these types of trends hold, it could mean that at least for the 
APSE supersonic engine studies a detailed stage-by-stage en- 
gine model may not be necessary. 


Future Plans 

In the future, the plan is as follows. For the engine, the plan 
is to further improve the engine model by incorporating more 
accurate information about its geometry and component per- 
formance maps as they become available. Develop a super- 
sonic inlet 1-D finite element model in MATLAB that can be 
integrated to the rest of the engine simulation to arrive at an 
end-to-end integrated propulsion system simulation. Run 
analysis to find out whether a component level lumped volume 
dynamic model is accurate enough or if a stage-by-stage com- 
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ponent model is needed. Assess the worst case disturbances 
for atmospheric wind gusts, pitch, yaw, and ASE modes exci- 
tation. Linearize the propulsion model and obtain TF for con- 
trol design, design controls (ref. 8), and integrate the controls 
into the nonlinear propulsion system simulation. Analyze 
thrust variations of the propulsion system due to expected 
upstream flow field disturbances. From these thrust variation 
results, the ASE team would be assessing couplings to the 
vehicle ASE modes and impacts to flow field disturbances, 
upstream of the propulsion system. Finally, integrate the pro- 
pulsion system with the vehicle ASE system model and refine 
the analyses of thrust variations, vehicle stability, and ride 
quality. As a separate effort, this relative high fidelity propul- 
sion system simulation will be utilized to design improved 
control laws in order to increase the efficiency of the propul- 
sion system. 

The Supersonics Project may later provide an engine for 
testing and validation of the models. In case the engine does 
not become available, or does not become available in time, 
the plan is to bound the problem. It is anticipated that the 
frequency responses of the propulsion system, in terms of the 
dynamic frequencies, will be primarily influenced by the par- 
ticular engine component geometries and those could possibly 
be estimated within an accuracy of few percent. The ampli- 
tudes of the responses would mainly be driven by the compo- 
nent performances, and these may also be estimated. Even 
though, the simulation results are rather sensitive to the accu- 
racy of the engine component maps. Similarly, the analysis 
could be bounded for the expected worst case upstream at- 
mospheric flow field disturbances and those flow disturbances 
that would be the result of excitations of vehicle ASE modes. 

Conclusion 

In this study a dynamic engine model is developed to en- 
able studies of thrust variations due to upstream flow field 
disturbances, such as wind gusts, pitch, and yaw, as well as 
disturbances due to excitation of the vehicle aero-servo-elastic 
modes. Ultimately, the goal is to integrate the propulsion sys- 
tem simulation with the vehicle structural dynamic model to 
create an integrated vehicle aero-propulso-servo-elastic model 
that will enable integrated control design studies and studies 


involving vehicle stability and ride quality. Another, goal in 
developing this control volume dynamic model is to enable 
controls design to improve propulsion system performance 
and efficiency. 

The engine component simulations are described in detail 
and pertinent results are given, including preliminary thrust 
variations due to the engine itself (i.e., without including the 
inlet at this point). Some comparative results are shown with 
previously published experimental data that validates the 
modeling approach. Enhancements to this engine model, in 
terms of more accurate geometry information and component 
performances are still in progress. The difficulty with these 
types of engine models is that more precise component per- 
formance maps and geometries are needed to achieve more 
accurate results and stable solutions. Future plans are also 
discussed that described the steps towards achieving an inte- 
grated aero-propulso-servo-elastic simulation to conduct the 
necessary research for supersonic vehicles. 
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